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Abstract. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction 
rates using only kinetic theory and fundamental molecular properties (i.e., no macroscopic reaction rate information) are 
extended to include reactions involving charged particles and electronic energy levels. The proposed extensions include 
ionization reactions, exothermic associative ionization reactions, endothermic and exothermic charge exchange reactions, 
and other exchange reactions involving ionized species. The extensions are shown to agree favorably with the measured 
Arrhenius rates for near-equilibrium conditions. 
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INTRODUCTION 

As we consider missions, both manned and unmanned, where the (re)entry velocities result in shock-layer tem- 
peratures on the order of tens-of-thousands of degrees Kelvin, the importance of electronic energy levels and reac- 
tions involving charged species becomes more pronounced. Although the treatment of these processes using the 
direct simulation Monte Carlo (DSMC) method is not a new development 1-3], most of these methods are based on 
measured equilibrium rates which are always questionable at high temperature, especially when applied to nonequi- 
librium conditions. The primary reason for the inadequate state of chemical reaction modeling is the difficulty in 
accurately measuring the internal energy state specific reaction rates to validate theoretical models in the tempera- 
ture range of interest. Therefore, most reaction rates are based on low-temperature equilibrium measurements and 
are fit to the most reliable measured data sets. In some cases, the uncertainty associated with the spread of these 
measurements exceeds one order of magnitude [4]. 

Molecular-Level chemistry models that predict equilibrium and non-equilibrium reaction rates using only kinetic 
theory and fundamental molecular properties (i.e., no macroscopic reaction-rate information) have recently been 
proposed[5-7] (in addition, a new vibrational relaxation model has been introduced that only relies on fundamental 
molecular properties and a single measure of vibrational collision number at a reference temperature as well as a 
new electronic energy relaxation model [8] following the kinetic-theory approach). Preliminary evaluations of the 
new chemistry model indicated that the resultant reaction rates are in very good agreement with measured Arrhenius 
rates for near-equilibrium conditions and with both measured rates and other theoretical models for far-from- 
equilibrium conditions[6]. The new chemistry model is extended in this paper to include various reactions involving 
charged species. 


DESCRIPTION OF NEW CHEMISTRY MODELS 

Bird[5], supported by Gallis[6], have proposed a set of molecular-level chemistry models based solely on fun- 
damental properties of the two colliding molecules: their total collision energy, their quantized vibrational energy 
levels, and their molecular dissociation energies. These models link chemical-reaction cross-sections to the energy- 
exchange process and the probability of transition between vibrational energy states. The Larsen-Borgnakke proce- 
dures and the principle of microscopic reversibility are then used to derive simple models for reactions. These mod- 
els do not require any macroscopic data and they function by seeking to balance the fluxes into and out of each state, 



thus satisfying microscopic reversibility. As mentioned in the introduction, preliminary evaluations [6] indicated that 
the resultant reaction rates are in very good agreement with measured Arrhenius rates for near-equilibrium condi- 
tions and with both measured rates and other theoretical models for far-from-equilibrium conditions. Since reactions 
involving only neutral species have been covered in References [5-7], these will not be discussed herein. A new 
methodology for accurately reproducing exothermic reactions has only recently been introduced by Bird[7], but the 
previous detailed balance procedure for exothermic reactions [5] appears to be sufficient if the molecules have simi- 
lar masses and the activation energy for the endothermic reaction is high enough (e.g. for real air). So, most reac- 
tions involving only atoms, molecules, and ions in air can be treated using the detailed balance procedure, but the 
reactions involving electrons will need to be revisited using the new methodology. 

Ionization 


The procedure for ionization is quite similar to that for dissociation. The difference is that ionization is linked to 
the electronic energy levels instead of vibrational energy levels as dissociation is. For ionization, the collision en- 
ergy is defined as 


E - E +E 

c trans,pair elec , particle 


(1) 


where the first term on the right is the relative translational energy of the colliding pair and the second term is the 
electronic energy of the particle being considered for ionization. If the collision energy is greater than the ionization 
energy of the particle, ionization occurs. 


Associative Ionization 

Associative ionization reactions require more care. These reactions can be thought of occurring as 

A + B -A A£* -> AB + + e (2) 

where AB* is an intermediate, excited molecule that may lead to the ionized species AB + plus an electron. The first 
step is similar to the recombination procedure described in Ref. [6]. If the AB* molecule is in the ground vibrational 
state after a Larsen-Borgnakke trial redistribution of the collision energy, the intermediate molecule is “formed.” In 
the second step, the dissociation energy is first added to the original collision energy. The new collision energy is 
then distributed among the internal energy modes of the AB* molecule and then this molecule is tested for ionization 
as described above. If the remaining collision energy after redistribution plus the electronic energy of the molecule 
is greater than the ionization energy of the molecule, the associative ionization reaction occurs. Since the reverse 
(exothermic) reaction occurs with an electron, the the new methodology of [7] should be used in place of the old 
detailed balance model of [5]. 

Endothermic and Exothermic Charge Exchange 

Charge exchange reactions are similar to the exchange reactions discussed in Ref. [6]. The endothermic reaction 

A + B + A + + B (3) 

takes place when the electronic energy level of the colliding neutral particle, A, after a trial Larsen-Borgnakke redis- 
tribution of the sum of the relative translational energy and electronic energy of the neutral particle is equal to the 
electronic energy level with energy just above the activation energy of the reaction. 

Since the exothermic charge exchange reaction occurs with atoms/molecules of similar mass, the detailed bal- 
ance model[5] is acceptable for use, but should be revisited using the new methodology [7]. The reaction takes place 
when the potentially formed neutral particle after a Larsen-Borgnakke trial redistribution of the total collision energy 
including the reaction energy is at the electronic energy level corresponding to that just above the activation energy 
of the endothermic reaction. 

Endothermic and Exothermic Exchange Reactions Involving Charged Species 

In exchange reactions involving charged species, the predominant mechanism is the breaking of a bond of one 
species (either neutral or charged) in order to exchange one atom to bond with the original atom, just as with an ex- 
change reaction involving only neutral particles. Therefore, the logic for performing this reaction involving charged 
species should be identical to that involving only neutral species. An endothermic exchange reaction takes place 



when the vibrational level of the colliding molecule after a trial Larsen-Borgnakke redistribution of energy resulting 
from the sum of the relative translational energy and the vibrational energy of the molecule is one level above the 
level corresponding to the activation energy E a 



(4) 


where i a is the first vibrational level above the energy threshold of the reaction, E a is the activation energy, fo is 
Boltzmann’s constant and O v is the characteristic vibrational temperature. 

Again, since the exothermic exchange reaction occurs with atoms/molecules of similar mass, the detailed balance 
model[5] is acceptable for use, but should be revisited using the new methodology [7]. An exothermic exchange 
reaction therefore takes place when the potentially formed molecule after a Larsen-Borgnakke trial redistribution of 
the total collision energy including the reaction energy is at the vibrational state corresponding to the activation en- 
ergy as in Equation (4). 


APPLICATION TO PARTICULAR REACTIONS 


The new chemistry models as discussed above are investigated by examining their ability to reproduce rates of 
chemical reactions for which measurements exist. The DSMC models are applied in a two-dimensional code where 
the results are accumulated over all cells to approximate a zero-dimensional simulation. The test gas is comprised of 
eleven species: O, N, O 2 , N 2 , NO, 0 + , N + , OE", N 2 + , NO + , and e. The simulations were run at varying temperatures 
between 10,000 K and 60,000 K at a number density of le23/m 3 in an adiabatic box 0.002 m on a side with 300,000 
particles. During the simulations, relaxation of the particles (including rotational, vibrational and electronic) was 
allowed to proceed as usual, but when a reaction was determined to take place, the number of reactions was ad- 
vanced by one but the simulators were left unchanged, so no energy was added to or taken away from the flow. In 
the figures included in this section, the lines are values from quoted rates in the literature while the symbols are the 
values that were sampled from DSMC. For a reaction A + B — > C + D, the forward reaction rate is calculated as 

dn. 

~r (5) 

k f = -&- 
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where the numerator is the change in number density of species A due to the reaction in the sampled time and the 
denominator is the product of the number densities of species A and B. 

The capability of these techniques to reproduce measured equilibrium reaction rates will now be investigated. 
However, the uncertainty of the currently available measured rates for atmospheric hypersonic flows can be large. 
Most of the experimental data on which reaction-rate sets are based were obtained for low temperatures (2000- 
7000K). Through careful study and comparison with numerical simulations, Park[4] suggests that the experimen- 
tally obtained reaction rates can be extrapolated up to 30,000 K. 

The uncertainty of these measured rates is not always reported and in some cases is unknown. The problem be- 
comes even more severe in the extrapolated portion of the tempera- 
ture range (T>1 0,000 K). Park[4] reports that some rates measured 
for the same reaction differ by more than one order of magnitude, 
which will be apparent in the figures that follow. 

Associative Ionization 

Due to their low threshold energies, associative ionization reac- 
tions play an important role in determining the level of ionization 
and structure of the flow field at hypersonic reentry conditions. 

There exists very little experimental or computational data in the 
literature regarding the associative ionization reactions in air. The 
associative ionization reaction of N + O has the lowest energy 
threshold, however the associative ionization of N + N and O + O 
are also very important mechanisms in the current application. The 
sampled DSMC reaction rates of these reactions are compared to 
measured rates from the literature[4, 9-11] in Figs. 1, 2, and 3, re- 
spectively. The discrete points represent the sampled DSMC rates, 
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FIGURE 1. Associative ionization rate of 
atomic nitrogen and atomic oxygen. 
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FIGURE 2. Associative ionization rate of atomic nitro- 
gen. 
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FIGURE 3. Associative ionization rate of atomic oxy- 
gen. 


while the lines represent the experimental measurements corresponding to the reference listed in the figure legend. 
As discussed earlier, there is a wide spread in the reaction rates reported in the literature. Even so, the sampled reac- 
tion rates for N + O and O + O (Figs. 1 and 3) go through the middle of the experimental scatter, while the sampled 
rates for N + N, while lower than the experimental data, are still within an order of magnitude of the lowest meas- 
ured rate. Given the assumed measurement uncertainty, the DSMC values are in good agreement with measured 
rates. 


Ionization 

Another important reaction mechanism for ionized, hypersonic reentry problems is the ionization of an atom or 
molecule by direct impact. There is very limited data in the literature for heavy particle impact ionization (such as N 
+ N — ► N + + N + e), so these reactions won’t be discussed in this work. There is, however, an appreciable amount 
of data related to electron impact ionization reactions, especially for atomic species. The electron impact ionization 
reaction rates of N and IV 2 are presented in Figs. 4 and 5, respectively, compared to values from the literature[4, 10, 
12-16]. As with the associative ionization data, there is a wide spread (almost three orders of magnitude) in the 
available reaction rate curves for the ionization of N, but the sampled DSMC data goes through the middle of the 
scatter. There weren’t nearly as many references for the electron impact ionization of IV 2 , but the comparison shown 
in Fig. 5 is still mostly within one order of magnitude (especially in the temperature range of interest for rarefied, 
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FIGURE 4. Electron impact ionization rate of atomic 
nitrogen. 
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FIGURE 5. Electron impact ionization rate of molecular 
nitrogen. 
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FIGURE 6. Endothermic charge exchange reaction N 2 + 
N + -► N 2 + +N. 
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FIGURE 7. Exothermic charge exchange reaction N 2 + + 
N^N 2 +N + . 


ionized, hypersonic reentry problems). Once again, the sampled DSMC values are in good agreement with the 
measured rates. 


Endothermic and Exothermic Charge Exchange 

In charge exchange reactions, a singly charged ion (such as N + or 0 + ) collides with a neutral atom or molecule 
and captures one of its electrons, thereby becoming a neutral atom. Only a small amount of energy is transferred to 
the electron donor, so the newly created neutral retains most of the original energy of the ion. These can therefore 
be an important mechanism for the creation of high energy neutral particles. The measured reaction rates from the 
endothermic and exothermic charge exchange reactions corresponding to N 2 + N + — > N 2 + + N are shown in Figs. 6 
and 7, respectively, compared to values from the literature [10, 11, 16, 17]. Again, the spread in the reported Ar- 
rhenius curves is quite large, but the sampled DSMC reaction rates are in good comparison with the reported values. 

Endothermic and Exothermic Exchange Reactions Involving Charged Species 

The remainder of relevant chemical reactions involve various exchange reactions involving charged species. As 
an example, the sampled endothermic and exothermic reaction rates related to O + NO + — > 0 2 + + N are presented in 
Figs. 8 and 9, respectively, compared to values from the literature[18]. The references for these types of reactions 
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FIGURE 8. Endothermic charge exchange reaction O + 
NO + -> 0 2 + + N. 
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FIGURE 9. Exothermic charge exchange reaction 0 2 + + 
N^O + NO + . 


are scarce, but for both the endothermic and exothermic reactions shown, the comparison between the sampled 
DSMC rate and the reported curve are within an order of magnitude. 


CONCLUDING REMARKS 

Extensions to the recently introduced kinetic-theory approach for computing chemical-reaction rates have been 
proposed for reactions including charged species. These reactions include ionization, associative ionization, endo- 
thermic charge exchange, and other exchange reactions involving charged species. These new models do not use 
measured macroscopic reaction rates to calibrate adjustable parameters. Instead, they make use of the principles of 
microscopic reversibility and molecular-level energy exchange to predict the probability that a chemical reaction 
occurs during a collision between to particles. 

For several reaction types, the DSMC models produce equilibrium reaction rates that are in good qualitative and 
quantitative agreement with the best available measured or extrapolated reaction rates. The uncertainty of the ex- 
periments do not allow an absolute assessment of the DSMC models through these comparisons. However, the dif- 
ferences between the most reliable reaction rates and the corresponding DSMC values are usually less than an order 
of magnitude, which is generally within the uncertainty of the measurements and often certainly within the scatter of 
equilibrium rates quoted in the literature. The advantage of the new kinetic-theory approach for computing reaction 
rates becomes apparent if the rate data is not available or if the application requires the extrapolation of the data to 
temperatures well outside the range of the measurements. 
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